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ABSTRACT 

Measurements of the neutral hydrogen gas content of a sample of 93 post-merger 
galaxies are presented, from a combination of matches to the ALFALFA.40 data release 
and new Arecibo observations. By imposing completeness thresholds identical to that 
of the ALFALFA survey, and by compiling a mass-, redsliift- and environment-matched 
control sample from the public ALFALFA.40 data release, we calculate gas fraction 
offsets (A fg a s) for the post-mergers, relative to the control sample. We find that 
the post-mergers have HI gas fractions that are consistent with undisturbed galaxies. 
However, due to the relative gas richness of the ALFALFA.40 sample, from which we 
draw our control sample, our measurements of gas fraction enhancements are likely to 
be conservative lower limits. Combined with comparable gas fraction measurements 
by Fertig et al. in a sample of galaxy pairs, who also determine gas fraction offsets 
consistent with zero, we conclude that there is no evidence for significant neutral 
gas consumption throughout the merger sequence. From a suite of 75 binary merger 
simulations we confirm that star formation is expected to decrease the post-merger gas 
fraction by only 0.06 dex, even several Gyr after the merger. Moreover, in addition 
to the lack of evidence for gas consumption from gas fraction offsets, the observed 
HI detection fraction in the complete sample of post-mergers is twice as high as the 
controls, which suggests that the post-merger gas fractions may actually be enhanced. 
We demonstrate that a gas fraction enhancement in post-mergers, relative to a stellar 
mass-matched control sample, would indeed be the natural result of merging randomly 
drawn pairs from a parent population which exhibits a declining gas fraction with 
increasing stellar mass. 
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1 INTRODUCTION 

Galaxies require gas for star formation. There exist several 
ways in which this gas can be acquired, such as feeding 
through intergalactic filaments, accretion of satellites and 
major mergers (e.g. Keres et al. 2005; Wang et al. 2011; 
L’Huillier, Combes & Semelin 2012; and see Sancisi et al. 
2008 and Sanchez Almeida et al. 2014 for reviews). There is 
compelling evidence that the gas reservoir must be continu¬ 
ally replenished if star formation is to continue (see Putman, 
Peek & Joung 2012 for a review and discussion of various 
refuelling mechanisms). For example, whilst the total stel¬ 
lar mass of galaxies on the red sequence increases signifi¬ 


cantly towards lower redshifts, both the mass in galaxies in 
the ‘blue cloud’ (Bell et al. 2004, 2007; Borch et al. 2006; 
Faber et al. 2007) and the mass density of HI in high column 
density QSO absorbers (e.g. Guimaraes et al. 2009; Noter- 
daeme et al. 2012; Zafar et al. 2013) exhibit no strong evolu¬ 
tion with redshift. This lack of evolution contrasts with the 
star formation rate (SFR) density, which evolves strongly 
over the epochs probed by these observations (e.g. Lilly et 
al. 1996; Madau et al. 1998; Hopkins & Beacom 2006). Re¬ 
fuelling of the galactic gas reservoir seems inevitable and 
is both widely predicted by simulations and supported by 
observations (e.g. Keres et al 2005; Oosterloo et al. 2007a; 
Sancisi et al. 2008; Brooks et al. 2009). 
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Galaxy mergers provide a natural way to deliver a fresh 
supply of gas to galaxies (although perhaps not sufficient to 
maintain current rates of star formation, e.g. Di Teodoro & 
Fraternali 2014). This delivery mechanism usually triggers 
new episodes of star formation in both the pre-coalescence 
phase (Barton, Geller & Kenyon 2000; Nikolic, Cullen & 
Alexander 2004; Alonso et al. 2006; Woods, Geller & Bar¬ 
ton 2006; Ellison et al. 2008a, 2013b; Li et al 2008; Scudder 
et al. 2012; Patton et al. 2013) and the final merger rem¬ 
nant (Wang et al. 2006; Karteltepe et al. 2010; Daddi et al. 
2010a; Genzel et al. 2010; Carpineti et al. 2012; Kaviraj et 
al. 2012; Shabala et al. 2012; Ellison et al. 2013a; Kaviraj 
2014). Such bursts of star formation may plausibly lead to 
a more rapid depletion of the gas reservoir. However, ob¬ 
servational results on this topic have yielded mixed results. 
Braine & Combes (1993) found the HI content of disturbed 
galaxies to be similar to that of unperturbed galaxies, and 
concluded that the HI content was unaffected by intereac- 
tions. In contrast, observations of HI in post-starburst galax¬ 
ies, which are implicated as merger-driven (Zabludoff et al. 
1996; Blake et al. 2004), find HI gas fractions that are low 
relative to late-types (presumably the merger progenitors) 
of the same stellar mass (Zwaan et al. 2013). However, com¬ 
pared to early types, the post-starbursts are relatively HI 
rich (Zwaan et al. 2013). 

In addition to supplying fresh gas, the merger may also 
re-distribute the gas, potentially making it unavailable for 
sustained star formation. For example, maps of HI in nearby 
galaxy pairs and post-mergers reveal long tidal gas tails and 
bridges that can extend up to a few hundred kpc (Hibbard & 
Yuri 1999; Wang et al. 2001; Koribalski & Dickey 2004). Gas- 
rich ellipticals, which may be merger remnants, sometimes 
exhibit large HI masses distributed in extended structures, 
but with surface densities that are insufficient to form new 
stars (Oosterloo et al. 2007b). Indeed, it has been suggested 
that the long lifetimes of such extended HI structures, which 
may have no (or very low surface brightness) optical coun¬ 
terparts (Hibbard & Yun 1999; Buyle et al. 2008) may be an 
effective way of identifying mergers long after optical distur¬ 
bances have faded (Hibbard & van Gorkom 1996; Holwerda 
et al. 2011; Lelli, Verheijen & Fraternali 2014). 

Despite these numerous individual studies, larger sam¬ 
ples, coupled with a comparison to isolated (i.e. non-merger) 
galaxies of similiar properties, are required to form a statisti¬ 
cal picture of the evolution of gas content during the merger 
sequence. In Fertig et al. (2015) we investigated the evolu¬ 
tion of HI gas fractions in the pre-merger phase of a sample 
of close galaxy pairs. Despite well documented SFR enhance¬ 
ments in this same sample (Scudder et al. 2012; Patton et 
al. 2013), Fertig et al. (2015) do not detect any decline in 
the HI gas fraction of close pairs. Indeed, in an earlier study 
based on a small sample of confusion-free HI measurements, 
Huchtmeier et al. (2008) suggested that mergers may be 
relatively Hl-rich. Casasola et al. (2004) present a similar 
study, but with over 1000 galaxies drawn from all interac¬ 
tion stages, from well separated pairs to coalesced remnants. 
HI gas fractions are again found to be enhanced by a fac¬ 
tor of ~ 10 in early type systems, but much more modestly 
over-abundant in late type galaxies. 

There are a few possible reasons why previous studies 
have not detected depletion in merger gas fractions. First, in 
order to obtain a complete view of gas depletion, it is impor¬ 


tant to consider both atomic and molecular gas, the latter 
of which may be expected to be more closely linked to star 
formation. One of the first studies of the combined atomic 
and molecular content of galaxy mergers was performed by 
Braine & Combes (1993), in which it was found that the 
molecular gas fraction is enhanced in interacting galaxies. 
Bettoni et al. (2001), from a study of the CO and HI con¬ 
tent of polar ring galaxies (which have likely formed through 
a past merger), also found gas fractions that are enhanced 
by an order of magnitude. Enhanced molecular gas frac¬ 
tions soon after a merger are predicted by simulations (Bckki 
2014), although eventually the gas is consumed by the trig¬ 
gered star formation. Stark et al. (2013) measured elevated 
molecular fractions in galaxies selected to have anomolously 
blue centres, which are frequently the signature of a recent 
merger. A similar result is found by Goncalves et al. (2014) 
who measure elevated CO fractions in their local sample of 
Lyman Break Analogs, many of which show merger proper¬ 
ties. Ultra-luminous infra-red galaxies (ULIRGs), which are 
often merger driven, also manifest very high CO gas frac¬ 
tions (Daddi et al. 2010b). Most recently, Ueda et al. (2014) 
have used a sample of 37 post-mergers to show that molec¬ 
ular gas disks are a common feature of merger remnants, 
indicating that these reservoirs have not yet been exhausted 
by merger driven star formation. 

A second potential reason for the lack of HI depletion in 
mergers is the relatively small consumption expected in the 
pre-merger (pairs) phase. Fertig et al. (2015) used a suite 
of binary merger simulations to show that, even at the clos¬ 
est separations, the HI gas mass is expected to be depleted 
by only ~ 0.04 dex. This is a relatively small effect that 
can be challenging to detect in small samples. Signatures of 
gas depletion may be more obvious in late-stage mergers, 
where simulations predict that SFRs will peak (e.g. Mihos 
& Hernquist 1996; Di Matteo et al 2007, 2008; Montuori et 
a. 2010; Torrey et al. 2012). In Ellison et al. (2013a) we con¬ 
firmed the prediction of high SFRs by combining a sample 
of galaxy pairs (with separations up to 80 kpc) and post¬ 
mergers. There is a steady increase in the SFR enhancements 
with decreasing separation, peaking in the post-mergers. We 
may therefore expect that the most rapid consumption of 
gas occurs after coalescence, depending on the rate of star 
formation and its duration. Indeed, simulations of galactic 
gas reservoirs of both atomic and molecular gas show an in¬ 
crease in gas fraction before a decline at later times (Bckki 
2014; Rafieferantsoa et al. 2014). The initial increase in HI 
gas fraction may be linked to the condensation of hot halo 
gas into the disk (Moster et al. 2011; Tonnesen & Cen 2012). 
High molecular gas fractions may follow from the increased 
central surface densities that precede starbursts in merg¬ 
ers. Indeed, several previous observational studies have mea¬ 
sured increased gas fractions in either atomic or molecular 
gas (Braine & Combes 1993; Bettoni et al. 2001; Huchtmeier 
et al. 2008; Stark et al. 2013; Goncalves et al. 2014). 

In this paper, we present a systematic study of HI 
masses in a moderately large sample of 93 recently merged 
galaxies. The sample is constructed from both new Arecibo 
observations, combined with publicly available HI masses 
measured as part of the Arecibo Legacy Fast ALFA 
(ALFALFA) project (Giovanelli 2005a,b; Saintonge 2007). 
Through a careful matching of the post-mergers with a sam¬ 
ple of isolated (non-merger) galaxies, it is possible to quan- 
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tify the enhancement or deficit in gas fraction associated 
with the merger process. Finally, we interpret our results 
with the aid of hydrodynamical merger simulations. 

We adopt a cosmology of Qa = 0.7, Qm = 0.3, Ho = 70 
km/s/Mpc. 


2 SAMPLES AND OBSERVATIONS 

2.1 The Ellison et al. (2013a) Galaxy Zoo 
post-merger sample 

Ellison et al. (2013a) identified and studied the optical prop¬ 
erties of 97 post-merger galaxies, selected from the SDSS, 
based on the visual classifications of the Galaxy Zoo project 
(Darg et al. 2010). In brief, the Galaxy Zoo visual identifi¬ 
cations are performed on some 900,000 galaxies in the SDSS 
DR6, as described in Lintott et al. (2008). For the compi¬ 
lation of their merger samples, Darg et al. (2010) addition¬ 
ally required 0.0005 < z < 0.1 and m r < 17.77. Further 
visual classification was performed by Ellison et al. (2013a) 
with the additional requirements that galaxies have an SDSS 
specclass=2 and that a Mendel et al. (2014) stellar mass 
must be available. The final sample is characterized by sin¬ 
gle late stage mergers that are highly disturbed (i.e. these are 
not galaxy pairs), with distinctive shells and tidal features 
and is therefore likely to be dominated by the remnants of 
gas-rich major mergers. Although it is not possible to know 
the exact mass ratio of the mergers that formed these rem¬ 
nants, simulations indicate that the mass ratios were likely 
closer than 1:4 (e.g. Lotz et al. 2010b; Ji, Peirani & Yi 2014). 


2.1.1 Matches with ALFALFA.40 

The Ellison et al. (2013a) post-merger sample was cross- 
correlated with the ALFALFA survey (Giovanelli et al. 
2005a,b). ALFALFA is a large survey conducted with the 
Arecibo telescope covering ~ 7000 square degrees of sky, 
with a frequency coverage between 1335 and 1435 MHz, i.e. 
with coverage of the 21 cm line extending to z ~ 0.06. The 
survey is conducted in a drift scan mode, such that sensi¬ 
tivity is naturally dependent on the galaxy distance. Gio¬ 
vanelli et al. (2005a) estimate the HI mass limit to be ~ 10 6 
M@ at a distance of 6.5 Mpc and 10' Mg at 20 Mpc. The 
40% data release (ALFALFA.40, Haynes et al. 2011) covers 
~ 2800 square degrees and contains ~ 15,000 extragalac- 
tic sources. For the cross-matching, we searched for AL¬ 
FALFA.40 detections within a 1.5 arcmin radius and a dif¬ 
ference in recessional velocity of less than 500 km/s for each 
post-merger. In practice, most matches were successfully 
made with offsets < 30 arcseconds. This procedure yielded 
8 matches, see Table [I] Of these 8, 7 are flagged as high 
quality detections in ALFALFA.40 (identified by the survey 
as code 1). In addition, there is one SDSS post-merger (ob- 
jID=587736808845082795) that is associated with an AL¬ 
FALFA.40 source that is flagged as a code 2 galaxjQ 


1 The ALFALFA catalog allocates a code 2 flag for galaxies with 
detections whose S/N<6.5. or with otherwise noisy spectra. In the 
case of objID=587736808845082795, the data cube is significantly 
affected by RFI, hence the classification as ‘code 2’, despite the 
relatively high formal S/N. 


For one object, objID=587736479206342888, the AL¬ 
FALFA.40 catalog reports a suspiciously low velocity width 
of 27 km s -1 . Upon inspection of the spectrum, the data 
were found to be very noisy with a possible second compo¬ 
nent in the emission. The original spectrum only measured 
one narrow peak resulting in a under estimate of the veloc¬ 
ity width. We applied further smoothing to the spectrum 
and re-measured the flux adding the second peak. From 
this spectrum we then re-computed the parameters for this 
galaxy with our own data reduction tools, and report the 
final adopted value in Table [T] 


2.1.2 Targeted Arecibo observations 

A further 37 post-mergers were observed during two runs 
(PI Fertig) with the Arecibo telescope on August 19 and 
December 28, 29 2013, see Table [2] The observations used 
the L-band wide receiver with the Interim Correlator and 
9-level sampling with 2 polarizations per board. As for the 
ALFALFA survey, we used the WAPP spectrometers with 
a total bandwidth of 100MHz centered on 1385 Mhz. The 
observations were designed to reach a 5<r H I mass limit of 
10% of the stellar mass. The minimum exposure time for the 
ON+OFF was set to be 4 minutes. For four of the sources the 
observing time was capped at 30 minutes, yielding a 5er HI 
mass limit of 50% of the stellar mass or better. For the data 
reduction, the spectra were Hanning smoothed and then re¬ 
duced using standard ON-OFF reduction. Five sources were 
observed on 19 August - two were detections and three were 
non-detections. An additional 32 sources were observed on 
28 and 29 December resulting in 10 more detections. Re¬ 
duced spectra of the 12 new detections are shown in Fig. 
[l] The final resolution ranges between 10 and 14 km s _1 , 
depending on smoothing. 

Table [2] lists various target and observational parame¬ 
ters, including exposure times, fluxes and the final derived 
HI masses for our Arecibo observations. HI masses (in units 
of solar masses) are calculated using the following formula: 


M{HI) = 2.356 x 10 5 x D 2 S 2 1 (1) 

where D is the distance to the source in Mpc ( cz/Ho , 
consistent with ALFALFA) and S 21 is the 21 cm line flux 
in Jy km s -1 . All detections have S/N > 5, with 11/37 of 
the spectra having S/N>6.5 (the requirement for a code 1 
classification in the ALFALFA.40 catalog). The limiting flux 
for non-detections is calculated from 


S 21 = S/N x 


RMS W 50 
w° s £ 1000 


( 2 ) 


where the RMS is measured from the spectrum and w S m 
is a smoothing width. This latter parameter is either W 50 / 2 O 
for W 50 < 400 km/s or 400/20=20 for W 50 > 400km/s. For 
non-detections, the full width at half maximum (FWHM) 
line width W 50 is calculated from the i-band Tully Fisher 
relation (Giovanelli et al. 1997). The 21 cm flux limits for 
the non-detections were then calculated from the inferred 
W 50 values inserted into equation [ 2 ] for a limiting S/N (sig¬ 
nificance) of 6.5 in order to be consistent with values quoted 
from ALFALFA. The non-detections in Table |T] therefore 
have the S/N fixed to a value of 6.5. 
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Table 1. HI masses from the ALFALFA.40 sample matched to post-mergers in the Ellison et al. (2013a) Galaxy Zoo sample. 


SDSS objid 

AGC # 

R.A. 

dec. 

redshift 

log M* 
(M 0 ) 

w 50 

(km s —x ) 

S21 

(Jy km s -1 ) 

RMS 

(mjy) 

S/N 

log M(HI) 
(M 0 ) 

588010880366477449 

212593 

172.62860 

5.8920 

0.03493 

10.26 

279i27 

1.86±0.12 

2.51 

9.90 

10.02i0.03 

587736479206342888 

714739 

228.22300 

10.1099 

0.03740 

8.90 

133i9 

0.83i0.13 

2.10 

8.60 

9.70i0.07 

587739132428157191 

727297 

237.06779 

25.5270 

0.04173 

10.06 

271i40 

1.19±0.10 

2.34 

6.90 

9.96i0.04 

587736543096799321 

253937 

226.32449 

8.1540 

0.03912 

10.01 

151±31 

0.96i0.07 

2.24 

7.80 

9.82i0.03 

587724233182412906 

113181 

28.62170 

14.1960 

0.05065 

9.64 

98i36 

1.35±0.06 

2.33 

12.90 

10.16i0.02 

587739382067822837 

727193 

234.00630 

25.5510 

0.03616 

10.44 

319±19 

1.33i0.10 

2.35 

7.10 

9.88i0.03 

587736808845082795 

245068 

223.05690 

12.0609 

0.05250 

10.48 

291i7 

2.01i0.09 

2.16 

12.20 

10.39i0.02 

588010880378404942 

8372 

199.99110 

5.8080 

0.02134 

10.51 

346il0 

2.89i0.10 

2.23 

15.50 

9.79i0.01 


The AGC number refers to the ALFALFA catalog number. W 50 is the FWHM of the 21 cm emission line, S 21 is the integrated 21 cm 
line flux. Derivations of the RMS and S/N are given in Haynes et al. (2011). Stellar masses are taken from the single Sersic fits of 
Mendel et al. (2014). For objID=587736479206342888 the spectral parameters (including W 50 , S 21 , the RMS and final derived HI 

mass) are determined from our own re-analysis of the ALFALFA data. 
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Figure 1. Spectra of the 21 cm emission lines for the 12 new detections in our post-merger sample (see Table[2]|. The vertical dashed 
lines indicate the velocity range over which the flux is summed and the vertical dot-dash line shows the redshift derived from the optical 
spectrum. The SDSS objID is given above each panel. 


Combined with the ALFALFA.40 matches, we therefore 
have a total of 45 HI mass measurements for the Ellison et 
al. (2013a) Galaxy Zoo post-merger sample, of which 20 are 
detections and 25 are upper limits. 


2.2 The Nair & Abraham (2010) supplementary 
sample. 

In order to expand our post-merger sample, we used the vi¬ 
sual classifications of ~ 14,000 galaxies from Nair & Abra¬ 
ham (2010). The Nair & Abraham classifications are also 
based on SDSS images, but from the DR4 and with slightly 
different magnitude and redshift thresholds from those of 
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Table 2. Target list and parameters for Aug. and Dec 2013 Arecibo observations. 


SDSS objid 

R.A. 

dec. 

redshift 

log M* 

(M 0 ) 

Time 

(ON+OFF) 

W 50 
(km s _1 ) 

S21 

(Jy km s -1 ) 

RMS 

(mJy) 

S/N 

log M(HI) 
(M 0 ) 

587731187807092815 

343.01901 

1.2500 

0.04854 

10.10 

1140 

190T10 

2.23T0.29 

2.33 

15.62 

10.35i0.06 

587731185669177359 

4.87780 

-0.6020 

0.03287 

9.32 

1800 

100±8 

0.57T0.07 

0.89 

14.14 

9.42T0.05 

588016840167653827 

117.37570 

18.9580 

0.04724 

10.37 

240 

351T16 

2.08T0.27 

1.82 

13.71 

10.31i0.06 

587728905026535460 

119.24500 

34.0929 

0.07073 

10.05 

1200 

342T12 

1.62T0.15 

1.02 

19.12 

10.54i0.04 

587738947741483063 

140.58290 

31.2859 

0.04871 

9.72 

1800 

227T15 

0.74T0.15 

1.20 

9.08 

9.88i0.09 

587738947205791935 

143.69679 

32.1430 

0.01486 

8.52 

1800 

122i6 

1.07T0.10 

0.95 

22.87 

8.99i0.04 

588848900971888657 

146.79910 

0.7030 

0.03055 

9.44 

1200 

170±9 

3.20T0.27 

2.07 

26.49 

10.11i0.04 

587739159266590725 

153.54409 

34.3429 

0.03758 

9.99 

240 

189T26 

1.08T0.15 

1.08 

16.23 

9.83i0.06 

587739305286303750 

172.32220 

35.5769 

0.03457 

10.96 

360 

369T15 

5.29T0.42 

2.66 

23.15 

10.43i0.03 

587741533323526200 

173.78129 

29.8910 

0.04624 

11.07 

240 

333T10 

1.25T0.16 

1.25 

12.24 

10.06i0.06 

587739507154288785 

185.06550 

33.6609 

0.02157 

10.46 

240 

414i6 

2.06T0.26 

1.98 

11.18 

9.62i0.05 

587734892216057938 

179.54049 

9.5360 

0.06985 

10.69 

240 

147T25 

0.71T0.12 

2.31 

5.66 

10.17i0.07 

588015510343712908 

20.46209 

0.8470 

0.09898 

11.15 

240 

438 

<0.91 

1.43 

6.5 

<10.59 

587730775492001995 

342.48558 

14.8750 

0.08993 

10.69 

540 

320 

<0.42 

0.81 

6.5 

<10.17 

587731185651810753 

325.23431 

-0.4780 

0.06816 

9.74 

1800 

185 

<0.51 

1.28 

6.5 

<10.01 

587738372206690585 

120.58889 

19.9939 

0.06771 

10.58 

240 

359 

<0.91 

1.65 

6.5 

<10.26 

587732470387703859 

129.57679 

33.5779 

0.06212 

11.18 

240 

470 

<1.45 

2.13 

6.5 

<10.38 

587732484342415393 

130.93750 

35.8279 

0.05393 

11.04 

240 

373 

<0.99 

1.75 

6.5 

<10.09 

587744729564512606 

133.78070 

15.2690 

0.07681 

10.99 

240 

523 

<1.07 

1.41 

6.5 

<10.44 

587738409785557168 

143.44760 

10.8120 

0.08558 

11.33 

240 

369 

<0.77 

1.38 

6.5 

<10.39 

587735349633351728 

148.30130 

13.1009 

0.07709 

11.34 

240 

374 

<0.83 

1.47 

6.5 

<10.33 

587738947747053602 

154.64019 

36.2239 

0.05412 

10.95 

240 

439 

<1.49 

2.33 

6.5 

<10.27 

587726032776265859 

159.13279 

2.3620 

0.05026 

10.30 

240 

279 

<1.15 

2.36 

6.5 

<10.10 

587734948595236905 

160.26559 

11.0959 

0.05297 

10.82 

240 

434 

<0.83 

1.32 

6.5 

<10.00 

587732702329045015 

160.29069 

6.8200 

0.03332 

10.57 

240 

344 

<0.88 

1.63 

6.5 

< 9.63 

588848900441112762 

160.79060 

0.3140 

0.09788 

11.01 

240 

590 

<1.77 

2.07 

6.5 

<10.87 

588017702386270301 

162.93029 

10.7570 

0.09734 

11.52 

240 

448 

<0.83 

1.28 

6.5 

<10.53 

588017703997276265 

163.71000 

12.0900 

0.08045 

10.55 

480 

244 

<0.53 

1.17 

6.5 

<10.17 

588848900981260449 

168.23460 

0.7260 

0.09687 

11.24 

240 

348 

<0.94 

1.73 

6.5 

<10.58 

587741533859414124 

171.14239 

30.0960 

0.05492 

10.74 

240 

257 

<0.71 

1.53 

6.5 

< 9.97 

587742013816897595 

175.26300 

22.0480 

0.06290 

9.91 

960 

297 

<0.65 

1.30 

6.5 

<10.05 

587742062680015040 

176.18339 

23.1620 

0.04839 

10.58 

240 

415 

<0.95 

1.57 

6.5 

< 9.98 

587742062680997918 

178.51669 

23.3400 

0.05093 

10.16 

240 

426 

<0.80 

1.29 

6.5 

< 9.95 

587726031175221368 

180.99819 

1.4110 

0.02230 

10.18 

240 

205 

<0.58 

1.39 

6.5 

< 9.10 

588017703469777068 

185.68820 

12.2939 

0.09209 

11.54 

240 

568 

<1.62 

1.96 

6.5 

<10.77 

588017730842853516 

201.19819 

8.2950 

0.05179 

10.94 

240 

293 

<0.85 

1.70 

6.5 

< 9.99 

587741603112157297 

201.27389 

27.5450 

0.07209 

10.85 

240 

405 

<0.94 

1.59 

6.5 

<10.32 


The AGC number refers to the ALFALFA catalog number. W 50 is the FWHM of the 21 cm emission line for detections and calculated 
from the i-band Tully Fisher relation for non-detections (Giovanelli et al. 1997). S 21 is the integrated 21 cm line flux. Derivations of the 
RMS and S/N are given in Haynes et al. (2011). The S/N is set to 6.5 for all non-detections. Stellar masses are taken from the single 

Sersic fits of Mendel et al. (2014). 


Galaxy Zoo: m g < 16 and 0.01 < z < 0.1. However, although 
the Nair & Abraham catalog is largely a subset of the much 
larger Galaxy Zoo sample, complementary classifications ex¬ 
ist, due to the different approaches taken by classifiers (the 
statistical output of non-expert crowd-sourcing vs. a single 
expert eye). In particular, the Galaxy Zoo classifications are 
biased towards the most highly disturbed mergers, which are 
obvious from a brief inspection of the SDSS postage stamp 
images provided to the public. Conversely, the Nair & Abra¬ 
ham ( 2010 ) catalog additionally identifies more subtle, and 
often low surface brightness, features such as faint shells and 
tidal tails, as we demonstrate below. 

We compile the supplementary sample of post-mergers 
in two steps. First, we make use of 411 potential late-stage 
post-mergers, galaxies with strong shells and tidal features 


and close companions with merger classifications not present 
in the original Nair & Abraham (2010) catalog. We supple¬ 
ment this with a query of the published catalog to identify 
galaxies with an interaction_flag > 2 , but with the pair_flag 
= 0, which yields 77 galaxies. These lists are not unique. 
Both lists are cross-matched with ALFALFA.40, combined, 
and then duplicates between the two, and with the Elli¬ 
son et al. (2013a) sample, are excluded. The final Nair & 
Abraham post-merger sample contains 48 galaxies with AL¬ 
FALFA matches, of which 23 have detections and 25 have 
upper limits (Table [3|. As for our targeted Arecibo sam¬ 
ple, we list the FWHM (W 50 ) of the 21 cm emission line 
for detections and the width derived from the i-band Tully- 
Fisher for non-detections, as described in the previous sec¬ 
tion. Derivations of the RMS and S/N are given in Haynes et 
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Table 3. HI masses from the ALFALFA.40 sample matched to post-mergers selected from Nair & Abraham (2010) 


SDSS objid 

AGC # 

R.A. 

dec. 

redshift 

log M* 
(M 0 ) 

W 50 

(km s -1 ) 

S21 

(Jy km s -1 ) 

RMS 

(mjy) 

S/N 

log M(HI) 
(M 0 ) 

587724232100282512 

100318 

8.77704 

14.3548 

0.03801 

10.02 

299i22 

2.23i0.10 

2.37 

12.1 

10.12i0.02 

587724232641937419 

110216 

20.011 

14.3618 

0.03104 

10.12 

83il3 

1.13i0.06 

1.94 

14.2 

9.66i0.02 

587724233720332316 


30.9965 

14.3104 

0.04269 

10.92 

384 

<1.28 

2.25 

6.5 

<10.00 

587726101750546619 

231606 

210.551 

4.58573 

0.04031 

11.02 

421i27 

1.68i0.12 

2.44 

7.3 

10.09i0.03 

587727221949530261 


357.996 

14.7511 

0.04668 

10.95 

408 

<1.33 

2.25 

6.5 

<10.10 

587727223020454084 


351.309 

15.2448 

0.04342 

10.40 

276 

<1.09 

2.25 

6.5 

< 9.95 

587728880873701381 


146.112 

4.49914 

0.04671 

10.58 

408 

<1.33 

2.25 

6.5 

<10.10 

587728881413390384 

201303 

152.655 

5.15067 

0.01367 

9.10 

134i9 

1.41F0.07 

2.32 

11.7 

9.11i0.02 

587728881418502205 

6049 

164.38 

5.69873 

0.05320 

11.16 

290i3 

1.98A0.11 

2.56 

10.1 

10.41i0.02 

587729158966345769 

222264 

194.588 

4.88565 

0.03612 

10.22 

355il0 

1.68i0.13 

2.67 

7.5 

10.00i0.03 

587729160042119249 


199.249 

5.65679 

0.03306 

10.77 

343 

<1.21 

2.25 

6.5 

< 9.76 

587730774424223763 


356.624 

14.8867 

0.05751 

10.77 

330 

<1.19 

2.25 

6.5 

<10.23 

587732053235794407 


118.396 

26.4624 

0.03645 

10.18 

190 

<0.90 

2.25 

6.5 

< 9.72 

587732577233469472 


154.101 

5.6954 

0.04899 

10.97 

325 

<1.18 

2.25 

6.5 

<10.09 

587732577774469219 

201673 

163.518 

6.72539 

0.02659 

10.38 

229i3 

1.96i0.08 

2.12 

13.6 

9.81i0.02 

587732578837201105 

191115 

138.302 

5.61087 

0.05846 

10.99 

206i7 

1.30i0.10 

2.71 

7.4 

10.30i0.03 

587732579385933949 

213056 

165.496 

8.10891 

0.03019 

9.76 

219il2 

1.78i0.09 

2.32 

11.6 

9.88i0.02 

587732702331535424 


166.047 

7.06904 

0.03255 

10.71 

274 

<1.08 

2.25 

6.5 

< 9.70 

587732702868865044 


167.004 

7.44229 

0.04111 

10.62 

374 

<1.27 

2.25 

6.5 

< 9.97 

587732703398527109 


150.393 

6.93203 

0.04892 

11.31 

428 

<1.40 

2.25 

6.5 

<10.16 

587732771053830198 


178.604 

9.60905 

0.03517 

10.55 

329 

<1.19 

2.25 

6.5 

< 9.80 

587732771053830214 

210905 

178.637 

9.54455 

0.03576 

9.70 

180ill 

1.06A0.09 

2.53 

7.0 

9.80i0.03 

587732772126785620 

210781 

176.792 

10.5177 

0.02137 

9.57 

100i9 

0.55i0.06 

2.21 

5.5 

9.08i0.04 

587732772130652231 


185.741 

10.5483 

0.02594 

10.54 

304 

<1.14 

2.25 

6.5 

< 9.52 

587734622163763258 

180008 

120.867 

25.1027 

0.02758 

10.08 

237i74 

0.96i0.09 

2.39 

5.8 

9.52i0.04 

587734891683119143 

220805 

188.681 

9.00471 

0.04306 

10.67 

286i65 

0.84i0.10 

2.44 

4.5 

9.85i0.05 

587734892216713241 


181.066 

9.3969 

0.03525 

10.64 

354 

<1.23 

2.25 

6.5 

< 9.82 

587734893287440431 

210497 

174.042 

10.0556 

0.02073 

10.41 

261i23 

2.87i0.11 

2.32 

17.1 

9.77i0.02 

587734894357446674 


165.312 

10.4811 

0.04222 

10.84 

346 

<1.22 

2.25 

6.5 

< 9.97 

587734948054696027 


151.83 

9.87247 

0.05735 

11.17 

452 

<1.48 

2.25 

6.5 

<10.32 

587735347489341442 

202805 

156.547 

12.5738 

0.03096 

10.14 

292il4 

1.77i0.10 

2.34 

9.9 

9.90i0.02 

587735348021624886 


145.944 

11.4564 

0.03906 

10.75 

270 

<1.08 

2.25 

6.5 

< 9.85 

587735348025884725 


155.81 

12.8733 

0.03295 

10.53 

262 

<1.06 

2.25 

6.5 

< 9.70 

587735349086912675 

181028 

126.597 

8.42377 

0.04594 

11.17 

205ill 

0.62i0.09 

2.67 

3.6 

9.77i0.06 

587736478124408901 


209.255 

11.6114 

0.03870 

10.72 

332 

<1.19 

2.25 

6.5 

< 9.89 

587736802935373840 


196.764 

13.637 

0.02727 

10.23 

312 

<1.16 

2.25 

6.5 

< 9.57 

587736813672202319 


237.478 

9.82864 

0.04587 

10.92 

404 

<1.32 

2.25 

6.5 

<10.08 

587738409251635215 

190666 

150.264 

11.4614 

0.03629 

10.34 

225±8 

3.17A0.09 

2.30 

20.5 

10.28i0.01 

587738410860806225 

190505 

146.827 

12.1585 

0.04751 

10.75 

216i6 

1.07i0.08 

2.25 

7.3 

10.04i0.03 

587738411401412648 


155.398 

13.7755 

0.04554 

10.91 

389 

<1.29 

2.25 

6.5 

<10.06 

588010359073931336 

201281 

151.836 

4.07931 

0.02859 

10.48 

358il0 

1.28i0.11 

2.48 

6.1 

9.69i0.03 

588010360148656136 

200217 

154 

4.95476 

0.03203 

10.02 

325i6 

1.72i0.10 

2.33 

9.2 

9.91i0.02 

588010879841992790 

231357 

201.04 

5.25972 

0.02489 

9.81 

229±4 

2.33i0.09 

2.31 

14.9 

9.83i0.02 

588010880371654746 

220306 

184.483 

5.86978 

0.02347 

9.72 

287i31 

2.01i0.11 

2.57 

10.3 

9.71i0.02 

588017569779744827 


203.101 

11.1064 

0.03146 

10.06 

302 

<1.14 

2.25 

6.5 

< 9.69 

588017702921633970 


159.469 

10.7583 

0.05251 

10.99 

432 

<1.42 

2.25 

6.5 

<10.23 

588017702940049506 

232401 

202.338 

11.3365 

0.04771 

10.40 

329il4 

1.20i0.12 

2.60 

5.7 

10.09i0.04 

588017703471611920 


189.9 

12.4389 

0.04086 

10.76 

383 

<1.28 

2.25 

6.5 

< 9.97 


The AGC number refers to the ALFALFA catalog number and exists only for detections. W 50 is the FWHM of the 21 cm emission line 
for detections and derived from the i-band Tully-Fisher for non-detections. S 21 is the integrated 21 cm line flux. Derivations of the 
RMS and S/N are given in Haynes et al. (2011) and are set to 2.25 mjy and 6.5 for non-detections (ALFALFA.40 limits). Stellar 

masses are taken from the single Sersic fits of Mendel et al. (2014). 

al. (2011) and are set to 2.25 mjy and 6.5 for non-detections 2.3 Control Sample 
in the ALFALFA.40 sample. 

In order to investigate differences between the post-mergers 
and the general population, it is necessary to define a con¬ 
trol sample of undisturbed galaxies (defined as having a 
Galaxy Zoo merger vote fraction of zero) with no close com- 
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Figure 2. The distribution of asymmetry (Ra) parameters in 
the 45 Galaxy Zoo post-mergers and 48 post-mergers from Nair &; 
Abraham (2010) with HI mass measurements, plus the combined 
control sample (red, blue dashed and black dotted histograms re¬ 
spectively) . Whereas the Galaxy Zoo post-merger sample has very 
few mildly asymmetric galaxies (Ra < 0.05), the Nair & Abra¬ 
ham sample includes galaxies with much more subtle asymme¬ 
tries. Both post-merger samples have significantly higher asym¬ 
metries than the control sample. 



Figure 3. Examples of highly asymmetric (Ra > 0.05) galaxies 
in the Galaxy Zoo sample (top row) and smoother (Ra < 0.05) 
galaxies from the Nair & Abraham supplementary sample in the 
lower row. The value of Ra is given at the top of each panel. 


panion (e.g. Ellison et al. 2013a,b). The control sample is 
compiled by matching non-merger galaxies simultaneously 
in stellar mass, redshift and environment Following our 
previous works, we compute local environmental densities, 


= < 3 > 

where d n is the projected distance in Mpc to the n th 
nearest neighbour within ±1000 km s . Normalized den¬ 


sities, S n , are computed relative to the median E n within 
a redshift slice ± 0.01. In this study we adopt n = 5. The 
tolerance for matching is 0.005 in redshift, 0.1 dex in stellar 
mass and 0.1 dex in normalized local density. If less than 
five matches are found, the tolerance is grown by a further 
A,s=0.005 in redshift, Alog M*=0.1 dex in stellar mass and 
A<5s=0.1 dex in normalized local density until the required 
number of matches is achieved. 

The control matching accounts for various intrinsic de¬ 
pendences and observational biases. Perhaps most impor¬ 
tantly, it has been shown that HI mass correlates well (al¬ 
though with significant scatter) with stellar mass (Fabello 
et al. 2011; Cortese et al. 2011; Catinella et al. 2012). Look¬ 
ing for deviations in the typical HI content therefore needs 
to take this scaling into account. Due to the fixed exposure 
time of the ALFALFA survey, the HI detection threshold 
is strongly distance dependent (Haynes et al. 2011), neces¬ 
sitating the redshift matching. Finally, the environmental 
matching mitigates biases due to contamination by other 
galaxies in the Arecibo beam (although manual rejection of 
cases with obvious beam contamination are performed as 
part of the catalog assembly), as well as any dependence of 
gas fraction on group or cluster membership, or other en¬ 
vironmental trends (e.g. Solanes et al. 2001; Kilborn et al. 
2009; Chung et al. 2009; Cortese et al. 2011; Fabello et al. 
2012; Catinella et al. 2013; Hess & Wilcots 2013). 


2.4 Sample comparison 

The final post-merger sample used in this paper is the com¬ 
bination of the Ellison et al. (2013a) Galaxy Zoo sample and 
the Nair & Abraham (2010) supplementary sample, and con¬ 
tains 93 post-mergers with HI masses (43 galaxies) or lim¬ 
its (50 galaxies). A summary of the number of targets in 
each sub-sample is given in Table [4] In Figure [2] we com¬ 
pare the asymmetries of the Galaxy Zoo and supplementary 
Nair & Abraham (2010) samples that have ALFALFA.40 
matches. The asymmetries are determined from the GIm2d 
bulge±disk decompositions of SDSS galaxies performed by 
Simard et al. (2011). In particular, we use the Ra parame¬ 
ter (measured in the r-band), which measures the fraction 
of light in the residual image (after subtraction of the bulge 
and disk) that remains after a 180 degree subtraction of it¬ 
self (see equation lib in Simard et al. 2002). Ra is therefore 
the residual asymmetric component. 

Figure [2] reveals that the Galaxy Zoo post-merger sam¬ 
ple is dominated by galaxies with strong asymmetries (Ra > 
0.05). Although some highly asymmetric galaxies do exist 
in the supplementary sampl^ it is clear that it contains 
a significant population of galaxies with Ra < 0.05. For 
comparison, the asymmetries of the control galaxies are also 
shown in Fig. [2] 60 per cent of the control galaxies have 
Ra < 0.02, a regime little populated by the post-mergers. 
Figure [3] shows some examples of galaxies with different 
asymmetries. In the top row we show five of the Galaxy 
Zoo sample with Ra > 0.05 and in the bottom row five of 


2 Unfortunately, without precise information on the progenitor 
morphologies, it is not possible to account for Hubble type in the 
matching process. 


3 Recall that the Nair & Abraham (2010) supplementary sample 
of post-mergers is checked against the Galaxy Zoo sample and 
duplicates excluded. The supplementary sample therefore has had 
some of its more asymmetrc galaxies excluded. 
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Figure 4. A comparison between the properties of the Galaxy Zoo post-merger sample and the Nair & Abraham supplementary sample 
(red and blue-dashed histograms respectively) for which we have HI measurements. Top left panel: Total stellar masses are taken from 
Mendel et al. (2014). Top right panel: Extinction and k-corrected Petrosian g — r colours. Bottom left panel: SFR enhancements. Bottom 
right panel: normalized local densities. 


Table 4. The number of post-mergers in each sub-sample 


Selected from: 

Ellison et al. 
(2013a) / Galaxy Zoo 

Nair & Abraham 
(2010) 

Total 

45 

48 

HI detections 

20 

23 

Above ALFALFA.40 

15 

16 

90% completeness 




the Nair & Abraham galaxies with Ra < 0.05. The value of 
Ra is given at the top of each panel. 

From simulations of galaxy mergers, it has been shown 
that the visual properties of post-mergers depend on 
timescale, the initial galactic properties, merger geometry 
and mass ratio (Lotz et al. 2008, 2010a,b; Ji et al 2014). 
For example, Lotz et al. (2010a) have demonstrated that 
higher gas fraction galaxies exhibit longer-lived tidal fea¬ 


tures and disturbed morphologies. Conversely, ‘dry merg¬ 
ers’, those with low gas fractions, will exhibit lower surface 
brightness, less dramatic, and shorter-lived features. Post¬ 
mergers with lower asymmetries (i.e. smoother morpholo¬ 
gies) may therefore be the result of lower gas fraction merg¬ 
ers, or they may simply be older. Trends of gas fraction with 
galaxy asymmetry are discussed later in this paper. 

The properties of the two post-merger samples are fur¬ 
ther explored in Fig. [4] Recall that some of the post-mergers 
in the Nair & Abraham (2010) catalog are also identified 
by Galaxy Zoo and duplicates are removed from the supple¬ 
mentary sample. The supplementary sample of post-mergers 
presented here represents only those galaxies uniquely iden¬ 
tified by Nair & Abraham (2010). The top left and right 
panels of Fig. [4] show that the stellar mass and Petrosian 
g — r colour distributions are similar between the two sam¬ 
ples. The bottom left panel of Fig. [4] compares, for galaxies 
classified as star-forming, the SFR enhancements, relative to 
th mass-, redshift and environment-matched control sample 
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described in Section [2.3| (A SFR, see Ellison et al. 2013a,b 
for more details). In brief, A SFR compares the SFR of the 
post-merger with that of its control sample: 
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A SFR = log SFRpm - log SFR contro i , (4) 

such that, for example, a factor of two enhancement in 
SFR yields a A SFR = 0.3. The median SFR enhancement 
is 0.41 for the Galaxy Zoo postmergers, but 0.62 for the 
Nair & Abraham sample. Moreover, there are no galaxies in 
the latter sample with a negative A SFR (lower SFRs than 
the matched control). However, the fraction of star-forming 
galaxies is significantly lower in the Nair & Abraham sample: 
35 per cent, compared with 52 per cent amongst the Galaxy 
Zoo post-mergers. 

There is also a difference in the local environments of 
the two samples, as shown in the lower right panel of Fig. [4} 
with post-mergers in the supplementary sample systemati¬ 
cally located in regions of lower galactic density. However, as 
noted above, the procedure for sample selection and remov¬ 
ing duplicates means that a comparison of galaxy properties 
between the two samples is not strictly based on physical 
properties (an analysis that we perform later in this paper). 
Nonetheless, the above statistics and Figs[2] [4] demonstrate 
that the two visually selected samples identify galaxies with 
different properties, despite being selected from the same 
imaging data. 


2.5 Completeness 

Of the 93 post-mergers presented in this study, 37 have new 
HI mass measurements from our Arecibo observations. It 
is interesting to compare the sensitivity of these new ob¬ 
servations to that of the ALFALFA survey. As described 
in Haynes et al. (2011), the completeness of the ALFALFA 
survey is basically a function of the expected FWHM of the 
21cm line, Wso- From equation 4 of Haynes et al. (2011), we 
have that the 90% completeness of code 1 sources (S/N>6.5) 
is given as 


log S 2 i 


f 0.5 log VFso - 1-14, if log Wso < 2.5 
1 log W50 — 2.39, if log W50 ^ 2.5 


(5) 


In Fig [5] we show the measured fluxes of ALFALFA.40 
code 1 sources, with the 90% completeness limit shown. 
We show in red the 45 measurements of Galaxy Zoo post¬ 
mergers (8 of which are taken from ALFALFA.40), and in 
blue we show the 48 measurements from ALFALFA.40 of 
the post-mergers from the Nair & Abraham (2010) cata¬ 
log. It is clear that some of our new observations are deeper 
than the completeness limit of ALFALFA, which is to be ex¬ 
pected from our observation strategy. In the following analy¬ 
sis, we will frequently consider only those post-mergers that 
lie above the completeness limit, in order to make fair sta¬ 
tistical comparisons between the gas content of the post¬ 
mergers and isolated (non-merger) galaxies. However, for 
the rest of this analysis, we will combine the supplementary 
Nair & Abraham sample with the Galaxy Zoo post-merger 
sample. We hereafter refer to this simply as ‘the post-merger’ 
sample; it contains 93 galaxies, of which 43 have HI detec¬ 
tions and 31 are above the ALFALFA.40 90% completeness 
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Figure 5. The 21 cm flux (521) as a function of the FWHM 
(W50). The black dashed line shows the 90% completeness of the 
ALFALFA.40 catalog, as given by Haynes et al. (2011), see Eqn 
[ 5 ] ALFALFA.40 (non-merger) galaxies are shown in grey. Post¬ 
mergers are shown in red and blue for the Galaxy Zoo (Ellison 
et al. 2013a) and Nair & Abraham (2010) samples, respectively. 
Of the 45 Galaxy Zoo post-mergers, new detections are shown 
with red triangles, new limits are shown with red arrows and 
matches to ALFALFA.40 are shown with red circles. All of the 
Nair & Abraham (2010) post-mergers have detections taken from 
ALFALFA.40, and are therefore all shown with circles. 


limit (Table [ 4 ]). Although we have shown that the properties 
of these samples differ, their properties overlap and the dis¬ 
tinction between them has no clean physical boundary. Later 
in this work, we will investigate how gas fractions depend 
on galaxy properties within this combined sample. 


3 RESULTS 
3.1 Gas fractions 

The HI gas fraction (hereafter simply ‘gas fraction’) is de¬ 
fined in this work as f gas = Mhi / M<L In Figure [6] we 
present the HI masses and gas fractions for the full post¬ 
merger sample. For comparison, we also show non-merger 
galaxies detected in the ALFALFA.40 sample. The isolated 
(non-merger) nature of these galaxies is determined by re¬ 
quiring that the Galaxy Zoo merger vote fraction be zero, 
and we further require that the ALFALFA ocode=I and 
Ngal=l, indicating that only one SDSS source has been 
matched. For clarity, only ALFALFA detections (not upper 
limits) are shown. Figure [6] shows that the galaxies in our 
post-merger sample are biased towards the high mass end 
of the galaxy distribution for which ALFALFA HI masses 
exist. Figure [6] also shows that some of our own Arecibo 


4 In theoretical works, it is more common to define f gas = Mhi 
/ M loin !: but we adopt the more frequently used definition in 
observational studies. 
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Figure 6. The HI mass (top panel) and HI gas fraction (bot¬ 
tom panel) as a function of stellar mass. Red points (and upper 
limits) represent measurements for the post-merger sample and 
small grey points are non-merger galaxies extracted from the AL¬ 
FALFA.40 data release (see text for details). 


observations are deeper than the typical ALFALFA.40 de¬ 
tections (see also Figure[5|, whose effective integration times 
are only 40 seconds (recall that ALFALFA upper limits are 
not shown to avoid crowding). 

In order to make a more fair comparison between the 
post-mergers and ALFALFA.40, Fig.[7]shows the same data 
as in Fig. [6] but with only those galaxies that pass the 
90% completeness threshold included. In order to test sta¬ 
tistically whether post-mergers exhibit distinct gas fractions 
from isolated galaxies, we will define a galactic ‘gas fraction 
offset’, A fgas, which quantifies a gas fraction relative to 
isolated (control) galaxies with otherwise similar properties. 
The A fgas metric is very similar in spirit (although opposite 
in sign) to the widely used HI deficiency (e.g. Haynes, Gio- 
vanelli & Chincarini 1984). The application of such ‘offset’ 
techniques has been previously used to effectively quantify 
differences in colour, SFR (e.g. Fig. |4|, metallicity and ac¬ 
tive galactic nuclei (AGN) frequency between galaxy pairs 
and non-mergers (e.g. Patton et al. 2011, 2013; Ellison et al. 
2010, 2011b; Scudder et al. 2012; Satyapal et al. 2014). 



Figure 7. As for Fig. [6] but only showing galaxies with HI 
detections that are above the ALFALFA.40 90% completeness 
limit. 


3.2 Gas fraction offsets 

The basic procedure for control matching is described in 
Section |2.3| although we now additionally require that both 
the post-merger, and galaxies in the control ‘pool’ lie above 
the ALFALFA completeness threshold, as defined in equa¬ 
tion [5] There are 31 post-mergers that fulfill this criterion 
(see Table [4| and can hence have robust gas fraction offsets 
quantified. Mergers that lie below the completeness thresh¬ 
old will appear to have spuriously low HI masses for their 
stellar mass, due to the lack of ALFALFA galaxies in this 
regime. As described in Section [2.3[ the control sample is as¬ 
sembled by simultaneously matching in stellar mass, redshift 
and local environment all ALFALFA galaxies that are also 
above the 90% completeness threshold. By definition, this 
is the majority, but not all, of the published ALFALFA.40 
sample. As described in Section [2.3| we require a minimum 
of 5 control galaxies for each post-merger; if this minimum is 
not reached, the matching tolerances are grown. In practice, 
there are typically 10 - 30 controls for each post-merger. 

From the matched control sample it is then possible to 
calculate a gas fraction offset defined as 
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Figure 8. Distribution of gas fraction offsets. In both panels, the 
solid black histogram shows the distribution of A f ga s for isolated 
ALFALFA galaxies and the red dashed histograms indicate the 
post-merger sample. In the top panel, the stellar masses of post¬ 
merger galaxies are taken from the Mendel et al. (2014) catalog. 
In the lower panel, stellar masses are taken from the MPA catalog 
(Kauffmann et al. 2003). 


A /'gas — log /gas,PM log /gas ,control 5 (6) 

where the control value is the median of the matched 
control galaxy gas fractions for a given post-merger. This 
calculation is only made for galaxies with detections above 
the ALFALFA.40 90 per cent completeness threshold. A 
value of Afgas = 0 therefore indicates that a galaxy has 
an HI gas content that is typical for its mass, redshift and 
local density. In contrast, positive and negative values of 
A f gas indicate gas-rich and gas-poor galaxies respectively. 
We can also calculate gas fraction offsets for isolated AL¬ 
FALFA galaxies, by again comparing their gas fractions with 
other isolated galaxies at the same mass, redshift and local 
density. 

The top panel of Figure [8] shows the distribution 
of Afg aa amongst isolated ALFALFA galaxies (black his¬ 
togram) and the post-merger sample (red dashed his¬ 
togram). For the post-merger galaxies that are above the 
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Figure 9. Distribution of mean gas fraction offsets as a function 
of projected separation for the pairs sample of Fertig et al. (2015, 
blue triangles) and the post-mergers (red circle located in the grey 
shaded box) presented in this paper. There is no indication for 
significant gas consumption at any stage during the merger. 


ALFALFA.40 completeness threshold the distribution of gas 
fraction offsets is very similar to the values observed in iso¬ 
lated galaxies. The mean gas fraction offset in the post¬ 
mergers is consistent with no merger-induced change in the 
atomic gas fraction: Af gas = 0.03 ± 0.03. 

Since the gas fraction combines both HI masses and stel¬ 
lar masses by definition, it is germane to consider whether 
the gas fraction results depend sensitively on the derivation 
of M*. We have adopted the single Sersic fits from Mendel 
et al. (2014), which include the numerous improvements to 
SDSS photometry detailed in Simard et al. (2011). How¬ 
ever, there are two possible concerns with these masses that 
are relevant to the study of post-mergers. First, a single 
Sersic fit does not account for flux in the asymmetric com¬ 
ponents of the light, which might be a significant fraction 
in post-mergers. Second, the Mendel masses do not include 
star formation histories with bursts. Although smooth star 
formation histories are a good average representation of the 
general galaxy population, this is not likely to be the case for 
mergers, which are known to exhibit bursts of star formation 
following close passages and at coalescence. For comparison, 
we therefore repeat the Af gas analysis with masses derived 
by Kauffmann et al. (2003) (hereafter ‘MPA masses’) that 
are based on total photometry and include bursty star for¬ 
mation histories. The comparison between values of A f ga s 
for the ALFALFA.40 sample and post-mergers derived from 
matching MPA masses, rather than Mendel masses, is shown 
in the lower panel of Fig. [8] Although the individual values 
of Afgas in the post-mergers are slightly different, the dis¬ 
tribution of values in the post-merger sample remains very 
similar to the isolated sample. The mean gas fraction off¬ 
set is also still consistent with the matched control sample: 
Af g as = -0.02 ± 0.04. 

Fertig et al. (2015) have recently used similar methods 
to ours to study the gas fraction offsets in galaxy pairs. 
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The only notable difference in their analysis (apart from the 
selection of pairs, rather than post-mergers) is that their 
control sample is composed of wide pairs, rather than truly 
isolated galaxies. The use of wide pairs as controls has been 
previously used in the literature (e.g. Scott & Kaviraj 2014) 
and Fertig et al. (2015) show that this is a robust approach 
for gas fractions. Therefore, although Fertig et al. (2015) use 
wide pairs, and this work uses isolated non-mergers, for the 
control sample, the methodology is otherwise identical and 
the results should be comparable. 

Fertig et al. (2015) found no difference between the gas 
fractions in the close and wide pairs, indicating that con¬ 
sumption is not significant in the pre-merger phase of the 
interaction. In Fig. [9] we combine the results of Fertig et 
al. (2015) with our post-mergers, in order to provide a com¬ 
plete view of HI gas consumption during the merger process. 
The blue points show the gas fraction offsets as a function of 
projected separation in the Fertig et al. (2015) pairs sample. 
The filled red point is the mean A f ga3 in the combined post¬ 
merger sample (Fig. [8|. None of the observed data points 
are located at A f gas < 0, indicating a lack of measured gas 
consumption throughout the merger stages probed by our 
samples. 

3.3 Detection fractions 

In order to make a fair comparison between the post-merger 
sample and our control sample (drawn from ALFALFA.40), 
it was necessary to impose a consistent detection threshold 
for both samples. Although this leads to a robust and fair 
comparison, the analysis of A f ga3 includes only one third of 
our sample galaxies. Therefore, whilst we can state that, for 
galaxies detected in HI above the ALFALFA threshold the 
typical HI gas fraction is consistent with isolated galaxies, 
we have not yet investigated the nature of the rest of the 
sample. It is possible, for example, that whilst some galaxies 
possess ‘normal’ gas fractions, the majority of post-mergers 
are actually depleted, and that this depletion is the main 
reason they are below our detection limit. 

In the absence of a deeper control sample, it is not pos¬ 
sible to quantify the depletion in either the detections that 
are below the ALFALFA.40 threshold (of which there are 
12), nor in the non-detections. However, we can compute 
the detection fractions of the post-mergers and compare it 
to that of the control, in order to judge whether the former 
are more frequently gas-poor. For this calculation, we re¬ 
make the control sample, identifying the single best match 
(minimized deviation in stellar mass, redshift and local den¬ 
sity) for each post-merger. This adjustment is necessary due 
to the different number of control galaxies per post-merger 
permitted by the original procedure. The detection fraction 
is defined as the fraction of galaxies above the ALFALFA de¬ 
tection threshold (equation [5j, such that this calculation is 
not affected by the higher sensitivity of some of our Arecibo 
observations. We find that the detection fraction of post¬ 
mergers is almost double that of the control (32±6 per cent 
versus 17±4 per cent). The result is not sensitive to the num¬ 
ber of controls required in the matching process (i.e. similar 
results are obtained with 5 matches). 

The high detection fractions measured amongst the full 
post-merger sample again support the lack of a dominant 
population of gas depleted galaxies. On the contrary, the 


factor of two enhancement in the post-merger HI detection 
fraction may indicate that either the post-merger galaxies 
have a higher tendency to contain some gas, or that their gas 
fractions have been preferentially boosted above the nominal 
ALFALFA completeness threshold. These possibilities can 
be reconciled with the gas fraction result (which yielded an 
average post-merger A f gas consistent with zero) by recall¬ 
ing that ALFALFA is a shallow survey that itself consists of 
relatively gas-rich members of the galaxy population. Our 
measurement of gas fraction offsets may therefore be an 
under-estimate of the true value^J Therefore, not only is 
there no evidence for gas consumption in the post-mergers 
based on direct comparison of the observed gas fractions in 
the control and post-merger sample (Fig|8j but our analysis 
of the HI mass detection fractions hint that post-mergers 
may actually be somewhat gas-rich. This possibility is ex¬ 
plored further below in Section [4] 

3.4 Dependence of Gas Fraction Offset on 
Asymmetry 

It is well known that post-mergers are most asymmetric soon 
after coalescence, and that tidal features will fade with time 
thereafter (e.g. Lotz et al. 2008), so that mergers selected 
from shallow imaging surveys will be biased towards recently 
coalesced galaxies (Ji et al. 2014). In addition to the effects of 
merger timing, the nature of the progenitor galaxies can af¬ 
fect the resulting asymmetries. For example, the simulations 
of Lotz et al. (2010a) show that smoother merger remnants 
result from more gas-poor mergers. 

In Fig. [2] we showed that there is a range of asymme¬ 
tries in our post-merger sample, with a bias towards the 
highest values in the Galaxy Zoo sample from Ellison et al. 
(2013a). Indeed, in Ellison et al. (2013a) we argued that 
these galaxies were only recently coalesced based on their 
low gas phase metallicities. It is therefore interesting to in¬ 
vestigate whether there is a trend of gas consumption with 
asymmetry, as might be expected if the least disturbed post¬ 
mergers are older. 

In the top panel of Fig. [10] we plot A f gas as a function 
of r-band asymmetry for the ALFALFA.40 complete post¬ 
merger sample. There is no obvious trend between asymme¬ 
try and gas fraction offset. In the lower panel we plot the dis¬ 
tribution of A fgas for high ( Ra ^ 0.07, red histogram) and 
low (Ra < 0.07, blue histogram) asymmetry post-mergers, 
relative to the median asymmetry of the sample. Based on 
the results of a Kolmogorov-Smirnov (KS) test, there is no 
statistical difference that the two asymmetry samples have 
different distributions of A f gas - 

An important caveat to the application of asymmetry 
metrics to mergers is their relative insensitivity to extended 
low surface brightness features. Adapting classic asymmetry 
indices to be more sensitive to such low surface brightness 

5 The same caveat is applied to measurements of other properties 
of mergers investigated with matched control samples, due to the 
intrinsic nature of detection thresholds. For example, Scudder et 
al. (2012) caution that the SFR enhancements in their sample 
of galaxy pairs may be under estimates if their control sample 
contains (isolated) galaxies which themselves have preferentially 
elevated SFRs that boost them across the detection threshold. 
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Figure 10. Top panel: Gas fraction offsets as a function of galaxy 
asymmetry for the post-merger sample. Bottom panel: Distribu¬ 
tion of A fgas for post-mergers split by asymmetry at R,\ = 0.07; 
the red and blue histograms show high and low asymmetry galax¬ 
ies respectively. 


features, particularly at large galactocentric radii, may help 
reveal a closer connection with changes manifest in galactic 
stellar and gaseous properties (e.g. Lelli et al. 2014; Pawlik 
et ah, in prep.). 


4 DISCUSSION 

One of the main goals of the work presented here is to tackle 
the question of whether the merger process significantly de¬ 
pletes the gas reservoir of the remnant. We have found no 
evidence of this in the measured HI gas fractions of our sam¬ 
ple of post-merger galaxies. Specifically, based on HI detec¬ 
tions in 31 galaxies for which we can identify a fair control 
sample from ALFALFA, we find that post-mergers have gas 
fractions that are indistinguishable from equivalent isolated 
galaxies. However, due to the shallow nature of the AL¬ 
FALFA.40 survey, the sample from which the controls are 
drawn is dominated by relatively gas-rich galaxies. Indeed, 
considering the full sample of 93 post-mergers, we find an HI 
detection fraction that is double that of the control sample of 
isolated galaxies. Therefore, contrary to finding evidence for 


gas consumption, our analysis of HI gas fraction enhance¬ 
ments suggests that post-mergers might be relatively gas- 
rich. A number of other observational studies have drawn 
similar conclusions (Braine & Combes 1993; Bettoni et al. 
2001; Huchtmeier et al. 2008; Stark et al. 2013; Goncalves et 
al. 2014). Previous simulations have proposed that mergers 
might renew their cold gas reservoirs through the cooling of 
hot halo gas (Moster et al. 2011; Tonnesen & Cen 2012). 
Indeed, recent cosmological simulations measured increases 
in the HI gas fraction of ~ 20 per cent within a few hundred 
Myr after a merger (Rafieferantsoa et al. 2014). 

However, despite the theoretical predictions for en¬ 
hanced gas fractions, there may be a more mundane way to 
create high gas fractions in recent post-mergers, originating 
from the declining gas fractions with increasing stellar mass, 
as shown in Figure[6] (and see also Catinella et al. 2012 and 
references therein). The simple addition of two randomly 
chosen galaxies in such a distribution will yield relatively 
high gas masses for a given stellar mass. We demonstrate this 
effect in Figure [TT] with a simple simulation. We draw 2000 
galaxies at random from ALFALFA.40, pairing them up by 
adding their stellar masses and gas masses to yield ‘measure¬ 
ments’ of fg as for 1000 fake post-mergers. For each of these 
fake post-mergers, we assemble a control sample matched in 
stellar mass and determine A f gas as described abov^] For 
reference, Figure 11 again shows (black histogram) the dis¬ 
tribution of A fgas for all ALFALFA.40 galaxies above the 
90% completeness limit that also appears in the two panels 
of Figure [8] The red dashed histogram shows the distribu¬ 
tion of the 1000 fake post-mergers, that are seen to be offset 
towards positive values with a mean A f gaB = 0.17 ± 0.01. 
This is a direct consequence of the combination of match¬ 
ing a control sample in stellar mass and the declining gas 
fractions with that same quantity. 

Interestingly, the histograms shown in Fig. [Tl] are quali¬ 
tatively similar to the results of Rafieferantsoa et al. (2014), 
who compare the gas fractions of recently merged galax¬ 
ies to non-mergers in a cosmological simulation and find a 
20 per cent increase in the former. We stress that the gas 
fraction enhancement seen due to this ‘summation’ effect 
depends on the gas fractions of the progenitor galaxies, and 
also on the distribution of gas fractions in their controls. If 
pre-merger galaxies have initially higher gas fractions, the 
resulting A f gas of the post-merger will also be shifted to 
more positive values. Our simulation of fake post-mergers 
is therefore illustrative of the qualitative offset towards en¬ 
hanced gas fractions that might be expected due to stellar 
mass matching. 

The elevated HI detection fraction (by a factor of two) 
in our post-mergers compared with isolated controls of the 
same stellar mass may plausibly be explained by this ef¬ 
fect. However, we see no net shift towards positive A f gas 
(Fig [8] as indicated might be expected by Fig. 11 As pro¬ 
posed above, the consistency between post-merger and con¬ 
trol gas fraction offsets may be due to a tendency for the 
controls themselves to be relatively gas-rich, which artifi¬ 
cially reduces the magnitude of the measured gas fraction 
enhancement. An alternative explanation for the lack of gas 


6 For the real post-mergers, we also match in redshift and envi¬ 
ronment, but this is not possible in this simulation. 
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Figure 11. Distribution of gas fraction offsets for isolated AL¬ 
FALFA galaxies (solid black histogram) and 1000 fake post¬ 
mergers (red dashed histogram), constructed through random 
draws from the ALFALFA.40 sample. Due to the declining gas 
fraction as a function of stellar mass, matching in stellar mass 
leads to a shift towards positive values of A f gas . 


fraction enhancement in the observations is that prompt 
consumption of the initial HI excess has already occured 
in the post-mergers due to high rates of star formation. In 
order to quantify the expected amount of gas consumption 
in the post-mergers phase, we turn now to a suite of binary 
merger simulations. 


4.1 Timescales for Gas Depletion 

Fig.[lI]demonstrates that, in a stellar mass matched sample, 
the trend of declining f gas with M* leads to a shift towards 
high gas fraction offsets when two galaxies are combined. 
Moreover, accretion of new cold gas is predicted by some 
merger simulations (e.g. Moster et al. 2011; Tonnesen & Cen 
2012), and also might be expected in generic scenarios of 
mixing a cool and ionized medium (Marinacci et al. 2010; 
Marasco, Fraternali & Binney 2012). However, Fig.[8]shows 
that there is, on average, no enhancement in the gas fractions 
in the observed post-mergers. In this section, we explore 
whether the contrast between Figs. [8] and 1 11 1 is actually an 
indication of gas consumption. 

A detailed assessment of gas consumption and recycling 
requires sophisticated hydrodynamical simulations that can 
adequately track the gas in its different phases (e.g. Dave et 
al. 2013), which is beyond the scope of this paper (but see 
Rafieferantsoa et al. 2014). However, we can make a rough 
assessment for a simplified case, based on the suite of merger 
simulations presented in Patton et al. (2013) using the model 
for radiative gas cooling, star formation, and feedback as 
summarized in Torrey et al. (2012). Such simulations will 
also provide an approximate timescale over which galaxies 
will return to their nominal scaling relations. 

The suite of binary mergers includes 75 orbital varia¬ 
tions of a merger with mass ratio = 1.4 x 10 10 


M@/ 5.7 x 10 9 Mg = 2.5, which represents the typical pair¬ 
ing in the SDSS sample of Patton et al. (2013). Two iso¬ 
lated galaxies, identical to each of the galaxies in the merger, 
which act as controls, are also simulated. By comparing the 
integrated SFR in mergers during the course of the simula¬ 
tion (which runs for a total of 6.5 Gyr), with the isolated con¬ 
trols, we can determine the typical additional stellar mass 
created in the merger (compared to the baseline stellar mass 
created in the controls), <5M. Although this simulation does 
not track the various phases of gas, we make the simplifying 
assumption that star formation is intially fed from gas that 
originated as HI. This assumption leads to a conservative 
upper limit on the gas consumption. 

In Fig. [12] we plot the accumulation of stellar mass 
throughout the 6.5 Gyr simulation for 3 out of 75 of our 
orbital variations, which include combinations of five ec¬ 
centricities and five impact parameters for 3 relative disk 
inclinations (Patton et al. 2013; Moreno et al. 2015). The 
red, green and blue lines show examples from the e, f and k 
suites, which have relative disk orientations of 60, 90 and 180 
degrees respectively (Robertson et al. 2006). The coloured 
dashed lines in the top panel shows the accumulated stellar 
mass in both galaxies in the merger, and the black line shows 
the mass build-up for identical galaxies evolved in isolation. 
The first (left-most) set of vertical dotted lines indicates the 
time of first peri-centre; this time is very similar for each of 
the 3 orbits. The second (right-most) set of vertical dotted 
lines shows the time of coalescence, which is more sensitive 
to the orbital parameters. The example e simulation (red 
line) is strongly affected by the first pericentric passage and 
significantly increases its stellar mass through triggered star 
formation after this time. However, it is little affected by 
coalescence. The opposite is true for the example k merger. 
Its anti-aligned disk means that it is little affected by the 
pericentric passage, but experiences a larger star burst at 
coalescence. The intermediate example, /, has modest star 
formation triggered at both peri-centre and coalescence (see 
Moreno et al. 2015 for a more detailed analysis of orbital 
and spatial dependences of triggered star formation). De¬ 
spite these very different histories, the total accumulated 
stellar mass is very similar for these three examples. This 
general tendency for small pericentre starbursts to be fol¬ 
lowed by large bursts at coalescence, and vice-versa, was 
also observed in the large suite of merger simulations run 


by Di Matteo et al. (2007). In the lower panel of Fig. 12 


we calculate the difference between the accumulated stellar 
mass in the mergers and that of the control galaxies. Again, 
the time at which the mass is formed is evidently depen¬ 
dent on the geometry, but converges to a similar value of 
~ 7 x 10 s Mg by the end of the simulation. Note that SM 
does not change much after 500 Myrs post-coalescence, in¬ 
dicating that even if there is a star burst at the time of the 
final merger, it is short-lived and gas consumption does not 
increase much further beyond the rate expected for the iso¬ 
lated galaxy. The range of SM (after 6.5 Gyr) amongst the 
full suite of 75 orbits spans 3 x 10 8 to 1 x 10 9 Mg, again, a 
relatively narrow range considering the different burst his¬ 
tories. 

Having established a typical value of SM = 7 x 10 8 
Mg, we assign HI masses (Mhi, 1 , Mhi, 2 ) to the galaxies by 
determining the best fit solution to the ALFALFA.40 data 
points shown in Fig. [7] 
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Figure 12. Top panel: Cumulative stellar mass formed in three 
example merger simulations from our suite of 75 orbital combi¬ 
nations of a 2.5:1 stellar mass ratio merger. The coloured dashed 
lines show the combined stellar masses made in both galaxies 
througout the 6.5 Gyr simulation. Red, green and blue lines show 
examples from the e, / and k orientations respectively. The first 
(left-most) set of vertical dotted lines indicates the time of first 
peri-centre. The second (right-most) set of vertical dotted lines 
shows the time of coalescence, which is sensitive to the orbital 
parameters. The black line shows the steady evolution of the 
isolated control galaxies. Bottom panel: The difference in stel¬ 
lar mass ( SM ) accumulated by the mergers relative to the control 
galaxies. 


log Mhi = 0.366 log M* + 6.231. (7) 

The gas fraction of the post-merger is therefore given as 

log fgas ,pm = [log(MHi,i + Mhi, 2 — 2 SM)] — 

[log(M* 5 i + M* )2 + 25M)]. (8) 

The control galaxy gas fraction f 9 as,c is determined from 

log fga,s,c = \og(M H I,c) - [log(M* 1 l + M*, 2 )]. (9) 

where Mhi,c is again determined from the best fit AL¬ 
FALFA.40 data for a stellar mass log(M*,i + M*, 2 ). In the 


case of no star formation (i.e. akin to the simulation of ran¬ 
dom pairings shown in the lower panel of Fig [8]) we find a 
mean A f gaB = 0.17. Inserting the appropriate numbers into 
equations [8] and [9] we find A f gas = 0.11 where star for¬ 
mation is depleting gas over the course of the simulation. 
Therefore, for the typical pair represented by our merger 
suite, the expected decrease in A f gas even several Gyr after 
coalescence is very modest, only 0.06 dex. For comparison, 
Fertig et al. (2015) used the same set of simulations to esti¬ 
mate an expected decrease in A f gas of 0.04 dex in the close 
pairs phase of the merger. 

Although based on simplified approximations, the re¬ 
sults of our simulations lead to similar conclusions to the 
more sophisticated treatment of gas fractions in mergers 
by Rafieferantsoa et al. (2014), namely that galaxy merg¬ 
ers may yield a medium-to-long term shift towards high gas 
fractions in their remnants, before a modest decline. The re¬ 
sults of these simulations indicate that the contrast between 
the expected shift to positive A f gas from the summation 


of galaxies (Fig. 11 1 and the lack of a shift in the observed 
post-mergers (Fig. 8| is not due to gas consumption. In¬ 
stead, we favour a scenario in which the matching of con¬ 
trol galaxies from the ALFALFA survey is biased towards 
gas-rich galaxies, which reduces the gas fraction enhance¬ 
ment. However, increased gas fractions (or gas presence) are 
nonetheless implied from the high HI detection fractions in 
the post-mergers. 


5 CONCLUSIONS AND CLOSING REMARKS 

HI mass measurements of 93 post-merger galaxies, visually 
selected from the Sloan Digital Sky Survey, are presented, 
including new Arecibo observations for 37 of the sample 
galaxies. The remaining HI masses are taken from the AL¬ 
FALFA.40 sample (Haynes et al. 2011). Through a carefully 
controlled matching procedure with isolated (non-merger) 
galaxies selected from the ALFALFA.40 data release, and 
taking into account differences in completeness thresholds, 
we calculate the gas fraction offset, A f gas , which quantifies 
the relative gas richness of each post-merger. We also calcu¬ 
late the gas fraction offset for all (non-merger) ALFALFA.40 
galaxies through the same matching procedure. Finally, we 
interpret our gas fraction results with the aid of a suite of 
75 binary merger simulations. 

Our conclusions are as follows: 

(i) The gas fractions in post-merger galaxies that are 
above the ALFALFA.40 completeness threshold are con¬ 
sistent with those measured in a mass-, redshift- and 
environment-matched control sample (Fig. [8|. Combined 
with comparable measurements in galaxy pairs (Fertig et 
al. 2015), we find that the HI gas fractions in interacting 
galaxies are consistent with their controls throughout the 
merger sequence (Fig. [9|, with no evidence for significant 
consumption of the neutral gas reservoir. 

(ii) There is no trend between the gas fraction enhance¬ 
ment and galaxy asymmetry (Fig. |10[ ). 

(iii) The HI detection fraction in the full post-merger 
sample is approximately double the detection fraction of the 
isolated, matched control galaxies. Therefore, in addition to 
finding a lack of gas consumption, our results suggest that 
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post-mergers may be, on average, more gas rich than isolated 
galaxies. 

(iv) Although high gas fractions in post-mergers have 
been predicted in several simulations (Moster et al. 2011; 
Tonnesen & Cen 2012; Bekki 2014; Rafieferantsoa et al. 
2014), we demonstrate that this enhancement is expected 
even in the absence of hot halo gas cooling, or the accretion 
of external gas. By drawing random pairs of galaxies from 
the isolated galaxies in the ALFALFA.40 sample we show 
that high gas fractions may arise naturally from the pair¬ 
ing of two galaxies drawn at random from the ALFALFA.40 
detections, due to the declining f gaa with increasing stellar- 
mass of the underlying population (Fig. 111. Indeed, the re¬ 
sults of this test indicate that for a stellar mass matched 
sample we might expect to see a shift towards positive A 
fgas in post-mergers. 

(v) Based on a suite of binary merger simulations, we 
compute a conservative upper limit for HI gas consumption 
for a galaxy interaction that is representative of the SDSS 
pairs sample (Fig 1 12[ ) . We find that, even several Gyrs after 
coalescence, the impact of star formation on gas consump¬ 
tion is small and that A f ga3 is expected to change very little, 
decreasing by only 0.06 dex (although the exact consump¬ 
tion is likely to depend on mass ratios, orbital geometries 
etc.). Gas consumption therefore does not appear to explain 
the lack of gas fraction enhancement in the observations. In¬ 
stead, we propose that our gas fraction offsets may be lower 
limits, due to the relatively gas-rich nature of the ALFALFA 
control sample. 


We conclude that there is no evidence for significant gas con¬ 
sumption in our post-merger observations; on the contrary, 
the post-mergers may actually be somewhat gas-rich. 


One remaining uncertainty in our investigation of the 
evolution of gas content in mergers concerns the molecular 
gas content. Measurement of the molecular component is not 
only necessary for obtaining a complete census of the galac¬ 
tic gas, but it is also critical for understanding the mode 
in which star formation is progressing. Highly star-forming 
galaxies at both high redshift and locally (e.g. the ULIRGs) 
are offset on the Kennicutt-Schmidt relation (Daddi et al. 
2010a; Genzel et al. 2010). Mergers have been predicted 
to also inhabit this offset sequence that exhibits ten times 
higher SFR at fixed gas (total) mass surface density (Renaud 
et al. 2014). 

Although numerous studies have investigated the de¬ 
tailed distribution of molecular gas in individual mergers 
(e.g. Yun & Hibbard 2001; Davis et al. 2013; Fernandez et 
al. 2014), a statistical study that combines HI and H 2 , with 
a well defined control sample is required in order to quantify 
general trends within mergers. One possible reason for the 
lack of HI depletion in our merger sample, despite elevated 
rates of star formation, could be that the stars form from the 
extant molecular gas reservoir. However, the previous stud¬ 
ies that have attempted to compare mergers and isolated 
galaxies have found that the molecular fraction is actually 
enhanced in mergers (Braine & Combes 1993; Bcttoni et 
al. 2001; Stark et al. 2013). Further hints on the impact 
of the molecular gas reservoir in mergers may be gleaned 
from observations of local UV luminous galaxies, or ‘Lyman 
break analogs’ (LBAs) selected to have properties similar 
to normal star-forming galaxies at high redshift (Heckman 


et al. 2005). Interestingly, the LBAs have many properties 
typical of local mergers, including high SFRs, low metallici- 
ties and large IR luminosities (Hoopes et al. 2007; Overzier 
et al. 2011). In a recent study of 6 LBAs, whose morpholo¬ 
gies and emission line kinematics again indicate interactions, 
Goncalves et al. (2014) found high CO fractions relative to 
‘normal’ low redshift galaxies (Saintonge et al. 2011). Whilst 
some enhancement in the CO fraction may be expected from 
the simple addition of reservoirs, as we found for the HI gas 
fractions, it is also expected that the strong gas flows in 
mergers, that ultimately drive star formation, will also in¬ 
crease the CO fraction (Braine & Combes 1993). Conversely, 
Sargent et al. (2014) conclude that starburst galaxies (which 
are not necessarily mergers) have CO fractions that are typ¬ 
ically 2-3 times lower than normal galaxies. 

The measurement of normal HI masses, but elevated 
molecular fractions in a small sample of mergers led Braine 
& Combes (1993) to propose that, contrary to depleting the 
gas supply, interactions might be responsible for the cooling 
of ionized halo gas, which leads to a central build-up of H 2 . 
Such a mechanism seems plausible theoretically, through the 
interaction of gas reservoirs (Marinacci et al. 2010; Moster 
et al. 2011; Tonnesen & Cen 2012; Marasco et al. 2012). Ob¬ 
taining CO measurements for the post-mergers in our sam¬ 
ple, and studying simultaneously the changes in the atomic 
and molecular gas components will therefore be an impor¬ 
tant complement to the study presented here. 
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